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Abstract 

The influence of the exit boundary conditions (vanishing first derivative of the velocity com- 
ponents and constant pressure) on the large eddy simulation of the fully developed turbulent 
channel flow has been investigated for equidistant and stretched grids at the channel exit. 

Results show that the chosen exit boundary conditions introduce some small disturbance 
which is mostly damped by the grid stretching- The difference between the fully developed 
turbulent channel flow obtained with LES with periodicity condition and the inlet and exit 
and the LES with fully developed flow at the inlet and the exit boundary condition is less 
than 10% for equidistant grids and less than 5% for the case grid stretching. The chosen 
boundary condition is of interest because it may be used in complex flows with backflow at 
exit. 
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Introduction 


Large eddy simulation of turbulent channel flows has been reported first by Kim and Moin 
[8] who considered a fully developed turbulent flow, for which periodicity conditions at the 
inlet and the exit of the computational domain can be assumed. The great advantages are 
that the problem of specification of the exit boundary condition can be avoided. The exit 
boundary condition for direct and large eddy simulation of turbulent flows plays a more 
important role with respect to damping than for laminar flows. However, in real flows the 
periodicity condition is more the exception than the rule. For example, for jet or impinging 
jet flows no periodicity condition can be used at the exit. Some research has been reported 
on the use of non-periodic exit conditions. 

Werner [18] used vanishing first or second derivative conditions for the velocity-components: 

d l U- 

M = O' , = 1 OT 2 TO 

With this von-Neumann-type-boundary-condition he calculated turbulent channel flows by 
using the Smagorinsky-Lilly [16] model. This kind of boundary condition is often used for 
laminar flows, where only spatially smooth fields of flow values like velocity, pressure with 
very small gradients exist. This means that further spatial development of flow is negligibly 
small. 

Werner [18] pointed out that a very smooth field can be represented in Fourier space as one 
wave with large wavelength and some waves with negligibly small wavelength. The difference 
form of vanishing first derivative corresponds exactly to waves with infinite wavelengths. The 
difference form of the vanishing second derivative with equidistant grids corresponds to a 
wave of infinite length and a wave with finite length, which is large in comparison to the 
grid distances. In contrast to laminar flows one can find very large gradients in turbulent 
flows, i.e., very rough fields. So there are also waves with short wavelengths, which are 
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reflected at the exit boundary. These reflections mean a disturbance of the results of the 
computational domain. Finally it should be mentioned, that vanishing higher derivatives 
as exit condition can lead to numerical instabilities [18]. An alternative class of boundary 
conditions are the so-called “N on- Reflecting- Boundary-Conditions” which were originally 
developed for hyperbolic equations [3-6]. Jin and*Braza [7] used such a condition for an 
elliptic problem. They calculated the transition of free shear layers using the two dimensional 
incompressible unsteady Navier-Stokes equations. 

The basic idea of their boundary condition formulation is the requirement that the waves 
of short length could pass through the boundary without being reflected. Basing on the 2D 
wave equation 

d 2 Ui 2 d 2 Ui 2 d 2 Ui 

~dW ~ Cx ~dx? ~ Cy ~dy* ~ 
they developed the following equation: 

dui dui d 2 Ui . 

-m- Ul -a^- v W = 0 

where u\ is the velocity in the main flow direction. 

Comparison of the coefficients of this equation with the Navier-Stokes equation provides the 
matched form of equation [3] . 

Another formulation of non-reflecting boundary condition is presented by Bestek [2] for the 
investigation of transition of boundary layers, from the solution of the velocity-vorticity 
formulation of Navier-Stokes equations. The second derivatives in the main flow direction 
(normal to the exit plane) are set to zero. So there is a damping which leads to relaminar- 
ization of the flow near the exit plane and reflection of disturbances is decreased. 
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Richter et al. [13, 14] have used “frozen turbulence” as the boundary condition in flows 
with pressure gradient. This technique is based on the assumption that the time variation of 
the turbulent fluctuations (the difference of instantaneous values of the flow variables from 
their avarage values) depend on the convective transport, i.e., the instantaneous values at 
the exit plane depend on the ones directly at the upstream. Richter et al. [13, 14] developed 
a relationship between the fluctuations at the exit plane and the plane directly upstream of 
the exit: 

^ ^0 |n = ^ (if, j', k) |n— An 

= Aa: 

’ j [At. Uc (j,k)} 

where u c is the characteristic convective velocity at the exit in x -direction and is set equal to 
< k ) >, where <> symbolizes the time averaged value. It is necessary to save enough 

time levels (An levels for getting | n -An) for implementation of this boundary condition. 
With the following equation 





§(i,j,k) = $[i,j,k) + $"{i,j,k) (6) 

the values at the exit plane can be defined. The average value $ is calculated by using a 
linear extrapolation of the variables from the interior points. 

The purpose of the present work is to perform LES for fully developed turbulent flows 
in a rectangular channel with an exit boundary condition and compare the results with 
computations with periodicity conditions in order to determine the influence of the boundary 
condition (see Fig. 1). The underlying idea is that for the channel with exit boundary 
condition the flow at the entrance is fully-developed turbulence. We would expect the 
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turbulence to remain fully developed and our chosen exit boundary condition should have 
negligible effect on the flow. 


periodic channel test channel 



Fig. 1: Schematic of the geometry 

As boundary condition we decided to use the vanishing first derivative condition like Werner 
[18] described above. One important aspect for this choice was the experience [9], that 
this condition in principle allows the calculation of complex flow structures like backflow 
at the exit plane caused by the entrainment. So it is interesting to investigate such a 
condition because there are no exact formulations for the LES of such kind of flows (like 
jets, impinging jets etc.) so far, where the flow is nonperiodic and a back flow appears at 
the nominal exit. Jin and Braza [7] showed indeed good results for the simulation of the 
2D Navier-Stokes equations for transitional free shear flow. But it is not known if their 
non- reflecting condition is suitable for calculating 3D turbulent flows by LES as well. Beside 
their boundary condition will not allow backflow at the exit plane. The same is true for the 
condition used by Bestek [2]. The “frozen turbulence” condition from Richter et al. [13, 14] 
is expensive because of the necessity of large computer time for sampling the variables at 
different time levels. 
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Numerical Method 


We consider a rectangular channel of length L, see Fig, 1. We use periodicity boundary 
conditions at x = 0 and x = L / 2 so that fully developed turbulent flow is simulated. Then 
the computation is further continued to x = L where the exit boundary condition is used. 
A comparison of the flow field at the first half and at second half of the channel would show 
the influence of the exit boundary condition (see Fig. 1). 


The basic equations are 


dui 

dxi 


0 


CO 


dui 

dt 


duiUj 

dxj 


duiUj 1 d 2 U{ 
dxj ' Re x dx\ 


For the subgrid scales the Smagorinsky- Lilly model is used (with C s = 0.1). 


UiUj — 2l/ijS{j 


( 8 ) 


Vi = (C s • A) 2 yj2SijSij 


( 9 ) 


4 - _ 1 dui duj 

2 dxj dxi 


( 10 ) 


Boundary Conditions 

No-slip conditions have been used on the walls for the convective terms. The no-slip condition 
for the diffusive terms were implemented by using the Schumann-assumption [15] in coupling 
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with the logarithmic law of the wall. Periodicity conditions have been used at x — 0 and 
x = i/2. At x ■== L we use 


duj 

dxi 


0 


( 11 ) 


and 


p = constant (12) 

Use of the coupled vanishing first derivative boundary condition with p = const, is based 
on an outflow boundary condition from Rannacher [12]. Our investigations showed that 
these conditions are suitable for the calculation of backflow caused by entrainment. The 
assumption of constant pressure is acceptable because at the exit plane the fluid leaves the 
computational domain towards the ambiant flow where a constant pressure, namely the 
ambient presssure, exists. The handling of backflow is obviously one of the great problems 
for simulation of impinging jet flow etc., see Benocci [1]. Some authors presented conditions 
where a so-called buffer layer near to the exit plane (Benocci, [1]) or the assumption of 
vanishing first derivatives in main flow direction for the base flow [10] is necessary. But this 
enlarges the computational domain and thereby the computational time. Because we want 
to investigate impinging jet flows by LES in future we decided to use Eqs. (11) and (12). 
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Method of Solution 


The basic equations have been solved by a fractional step method following Kim and Moin [8] 
using Adams-Bashforth and Crank-Nicholson difference schemes and a SIP solver of Stone 
[17] for the numerical solution of the pressure Poisson equation. 

Results and Discussion 

The Reynolds number based on the half of the channel height is 

Re _ = 2800, (Re T = 180) (13) 

u m is the average velocity in main flow (a;—) direction. Each half of the channel was dis- 
cretized with 34x34x18=20808 grids. The nondimensional length in the ^-direction was 8.5 
for each half, in the y-direction 4 and in the ^-direction 2 (Aa;=0.25, A y = Az = 0.125). 

In addition a modification of the boundary condition was done. 


periodic channel test channel 



Fig. 2: Schematic of the computational scheme 



The length of the second half of the channel, i.e., the part where we wanted to investigate 
the influence of the boundary condition, was increased to L = 9.145 by stretching the last 
8 grid cells in ^-direction steadily by 10% of the grid distance of the upsteam neighbouring 
cell (see Fig. 2). With this stretching, the numerical viscosity is increased. This stretching 
is probably equivalent to a buffer layer. The computations will also show if and how the 
nonequidistant grids would influence the results in comparison to the case of equidistant 
grids. 

As results we show first the distributions of different root-mean-square(rms)-values. For the 
investigation of the influence of our boundary condition we divided the second half of the 
channel into two parts (each part with 17 grid cells), see Fig. 2, so that one can observe if 
and how the influence of the exit plane is transmitted in the upstream direction. Figures 3, 4 
and 5 compare the distribution of the different rms values in the three directions. Obviously 
the rms values for the channel with our boundary conditions are higher than that of the 
periodic channel. The graphs of the channel half with our boundary conditions are nearly 
identical. This indicates an influence of the exit boundary condition through the whole 
channel. The graphs of the test cases show (in the periodic case) nearly constant deviation 
in the channel core and are identical in the wall region. In Table 1, the local maxima (max) 
and minima (min) of each rms-component of the periodic and the two parts of the test 
channel are compared. The differences of the each part of the second channel expressed in 
percent are related to the periodic channel. The values of the two parts of the test channel 
differ at both maxima and minima. But the largest difference, the rms- value in y-direction, 
is only 10.27%. All the other differences are lower. Our boundary condition introduces a 
disturbance in the whole test channel, which is expressed by higher rms values. But the 
effect of this disturbance is relatively low. Both parts of the test channel are influenced 
equally. 
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It should be mentioned, that the differences of the maximum rms values are smaller than 
the differences of the minimum values, see Table 1. This means that the locations of the 
flow where the fluctuation level is high are less influenced than locations where the level is 
low. 

Figures 6-8 show the distribution of rms values for the case of grid spreading. The graphs of 
the first part of the test channel are above the ones of the periodic channel. The values of 
the second part are equal or, even lower at the maximum value than the ones of the periodic 
channel. This indicates a decrease of the fluctuations because of the damping effect of the 
grid spreading. Table 2 shows the maximum and minimum rms values and the differences 
of the two parts of the test channel from the periodic solutions. The mentioned tendencies 
can be confirmed again. The rms values of the first part are higher than that of the periodic 
channel, but lower than the ones for the case of equidistant grids, (see in Table 1). All 
the values of the second part are even lower than the values of the periodic channel. This 
means, that in case of grid spreading there is also an influence of the boundary condition. 
But compared to the case of equidistant grid a decrease of the differences with respect to the 
periodic calculation could be achieved. The reason is the increase of the damping effect by 
numerical viscosity because of the grid spreading. This spreading influences very strongly 
the calculation of the rms values of the second test channel half. In this part, the damping 
effect is larger. This explains the low rms values here. Besides, because of this damping the 
disturbance caused by the boundary condition cannot be reproduced as in case of equidistant 
grids. There is a decrease of the values in the first test channel part. 

In Figs. 7 and 8 the Reynolds shear stresses at x — z plane (see Figures 1 and 2) are shown. 
Table 1 compares the maximum and minimum values of the Reynolds stresses. One can 
point out that the simulation using the present boundary condition give acceptable results. 
In case of equidistant grids the maximum difference from the periodic case for the local 


9 



maxima and minima is -6.17% (see Table 1, the negative sign results from the definition of 
the direction of the axis-system) . The comparison of the averaged centerline velocity confirm 
this conclusion, see Table 4. 

In case of grid spreading the results are even better, see Fig. 8 and Table 4. Here the 
maximum difference for the comparison of the local maxima and minima is only 2.40%. The 
comparison of the centerline velocities also show very good results, see Table 4. 

Conclusion 

The influence of coupled exit boundary condition of a vanishing first derivative on LES 
of turbulent channel flows for the velocities and constant pressure was investigated. A 
comparison of the results with the case of using periodic boundary condition was done. 
Besides, a modification - continously stretching of 8 grid cells near the exit plane - of the 
selected boundary condition was presented. 

The results confirm the suspicion, that the disturbance caused by the boundary condition at 
the exit plane is transported in the upstream direction. The influence of the chosen boundary 
condition leads to slightly increased turbulent fluctuations. This behavior is expressed by 
higher rms values. 

The results with grid spreading, because of the increased damping of non-equidistant grids, 
show better agreement with fully developed results. 

It was also found that the disturbance caused by these boundary conditions are small where 
the flow fluctuates strongly. So one can expect that this boundary condition should give 
good results for other kinds of flow especially with high levels of turbulent fluctuations like 
impinging jets. 
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v-RMS era u-RMS 


• periodic channel 
a 1 . part of the second channel 

□ 2. part of the second channel 



. 3: u-rms, equdistant grid 


• periodic channel 
a 1 . part of the second channel 

□ 2. part of the second channel 



Fig. 4: v-rms, equdistant grid 
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u-RMS (K) w-RMS 


• periodic channel 
a 1 . part of the second channel 

a 2. part of the second channel 



. 5: w-rms, equdistant grid 


• periodic channel 
a 1 . part of the second channel 

□ 2. part of the second channel 



Fig. 6: u-rms, non-equdistant grid 
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w-RMS (Jq v-RMS 


▲ 

□ 


periodic channel 

1 . part of the second channel 

2. part of the second channel 



.7: v-rms, non-equdistant grid 


• periodic channel 
a 1 . part of the second channel 

o 2. part of the second channel 



Fig. 8: w-rms, non-equdistant grid 
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• periodic channel 
a 1 . part of the second channel 
•0 2. part of the second channel 



Fig. 9: Reynolds stresses (uw-component in xz- plane) 


• periodic channel 
a 1 , part of the second channel 

□ 2. part of the second channel 



Fig. 10: Reynolds stresses (uw-component in xz-plane) 
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■ 

periodic 

channel 

1. part 
testchannel 

2. part 
testchannel 

difference 
of 1. part 

difference 
of 2. part 

u-rms 

min 

1.028 

1.075 

1.067 

4.85% 

3.88% 

max 

2.112 

2.146 

2.154 

1.61% 

1.99% 

v-rms 

min 

0.584 

0.644 

0.643 

10.27% 

10.10% 

max 

0.883 

0.924 

0.928 

4.64% 

5.10% 

w-rms 


0.549 

0.604 

0.601 

10.02% 

9.47% 


0.635 

0.666 

0.672 

4.88% 

5.83% 

Reynoldsstress 

) 

mm 

-0.519 

-0.543 

-0.551 


-6.17% 

max 

0.529 

0.554 

0.550 

g 

3.97% 


Table 1: Comparison of local maxima (max) and minima (min) for different lD-ensemble- 
averaged values, equidistant grid (averaging was done in x- and y-direction; the differences 
of each part of the second (test-) channel expressed in percent are refered to the values of 
the periodic channel) 




periodic 

channel 

1. part 
testchannel 

2. part 
testchannel 

difference 
of 1. part 

difference 
of 2. part 

u-rms 

m 

0.992 

1.016 

0.989 


-0.30% 


2.117 

2.139 

2.081 


-1.70% 

v-rms 


0.587 

0.620 

0.584 

g^^gggi 

-0.51% 



0.921 

0.859 


-3.59% 

w-rms 

B 

J||||g 

0.578 

0.542 

6.06% 

-0.55% 


| 

0.671 

0.629 

2.76% 


Reynoldsstress 

min 


-0.540 

-0.526 

-0.93% 

1.68% 

max 

0.513 

0.518 

0.501 

0.97% 

-2.40% 


Table 2: Comparison of local maxima (max) and minima (min) for different lD-ensemble- 
averaged values, non-equidistant grid (averaging was done in a:- and y-direction; the differ- 
ences of each part of the second (test-) channel expressed in percent are refered to the values 
of the periodic channel) 
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periodic 

channel 

1. part 
testchannel 

2. part 
testchannel 

difference 
of 1. part 

difference 
of 2. part 

Centerline- 

velocity 

18.65 

' 

18.57 

18.54 

-0.43% 

-0.59% 


Table 3: Comparison of the averaged Centerline- velocities, equidistant grid 



periodic 

channel 

1. part 
testchannel 

2. part 
testchannel 

difference 
of 1. part 

difference 
of 2. part 

Centerline- 

velocity 

18.58 

18.51 

18.48 

-0.53% 

-0.53% 


Table 4: Comparison of the averaged Centerline-velocities, non-equidistant grid 
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